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Abstract 



1 Introduction 

X 

' Recently, there has been a great deal of interest in Ginsparg- Wilson fermions |^ , because it is now realised 

that they allow the calculation of chirally symmetric physics without having any doubling problem ^, ^, Q . 
We also know that because of its chiral properties the Ginsparg- Wilson fermion matrix is automatically 
an 0{a) improved action in exactly the same sense that the clover action is an improved action. 

If we are interested in going beyond spectrum calculations to compute improved matrix elements, 
for example for structure functions and decay constants, we also need to know how to improve fermion 
operators. When we compute forward hadronic matrix elements, it is enough if the operators are improved 
for on-shell quantities. However, some methods of doing non-perturbative renormalisation ^ require 
calculations of off-shell Green's functions, with a virtuality large enough that we can reasonably compare 
with continuum perturbation theory. To do this well, we would like to be able to remove 0{a) effects 
from off-shell Green's functions too. 

In this paper we find several ways of improving operators for the Ginsparg- Wilson action, both for 
on-shell and off-shell Green's functions. The paper is organised as follows. In sect. 2 we first study 
the fermion propagator, and in sect. 3 we test the results in the free case for a particular realisation of 
Ginsparg- Wilson fermions. In sect. 4 we then turn to our main subject, operator improvement. The 
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major outcome is that, in the simplest form, the improved operators do not require any coefficients to 
be tuned non-perturbatively. In sect. 5 we briefly comment on the Ward identities, and in sect. 6 we 
conclude. 



2 Fermion propagator 

In this section we review some known results concerning the Ginsparg- Wilson propagator, as preparation 
for our consideration of bilinear fermionic operators. The basic Ginsparg- Wilson condition is |l| 

DgwIz+I^Dgw ^ aDcwl^Daw, (1) 



where Dqw is the Ginsparg- Wilson fermion matrix. The fermion matrix does not anti-commute with 75, 
but there is nevertheless a form of chiral symmetry present Q. 

Let us now look more closely at the Ginsparg- Wilson matrix. From the matrix Dqw we can define a 
related matrix 0] 

KGw=(}~'^DGwy^ Dgw. (2) 
We will discuss later what happens when Dqw has an eigenvalue equal to 2/a. This implies 

Daw = (1 + '^Kgw) ' Kgw- (3) 

The eigenvalues of Dqw lie on a circle of radius 1/a and centre 1/a, while the eigenvalues of Kqw lie on 
the imaginary axis. The relationship between the eigenvalues of the two matrices is shown in Fig. |l|. If 
we substitute eq. (^) into the condition (Q) we find that 

Kgw 75 + 75 Kgw = 0. (4) 



The chiral properties of Kgw are even closer to those of the continuum Dirac operator than Dgw- 
The fermion propagator we really want to use is the propagator calculated with Kgw- The propagator 
calculated from Dgw has contact terms of 0(a), which mean that if we Fourier transform it we will 
find 0(a) lattice artifacts when we are off-shell. A propagator calculated from Kgw contains no contact 
terms: it satisfies chirality even at zero distance. Therefore it will be improved off-shell too. 

We want to write down improved Green's functions for the massive as well as for the massless theory. 
When we consider massive fermions, we will use the fermion matrix |^] 

M = ^1 - ^amo^ Dgw + m-o • (5) 

This is the usual choice when working with Ginsparg- Wilson fermions, and it has the advantage of giving 
a simple linear relationship between bare and renormalised fermion masses. Most other ways of adding 
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Figure 1: The spectrum of the Ginsparg- Wilson matrix Dqw (the circle) and of the associated matrix 
Kgw (dotted vertical line). The dashed line links an eigenvalue of Dqw to the corresponding eigenvalue 
of Kgw- 

a quark mass lead to non-linear relationships of the form mji oc mo(l + arriQbjn + •••). However, despite 
the slightly more complicated form of the answer, the improvement schemes we suggest in this paper also 
work with minor modifications for other choices of M ^ for example M — Dgw + ttlq. 

The fermion propagator 



still contains contact terms violating chiral symmetry.]^ The fermion propagator with the correct chiral 
symmetry properties is the propagator corresponding to Kgw, i-e. 



We will always use ★ to denote improved quantities. The improved Green's functions have the correct 
chiral properties, and we therefore know that they are free of 0{a) discretisation errors. This is because 
any 0{a) term in a Green's function has the opposite chirality to the leading term, so checking chirality 
is a simple way of testing for 0{a) discretisation errors. The improved propagator (^ can be calculated 

^ We will be considering Green's functions involving charged or coloured fermions. These are, of course, only defined 
if the gauge is fixed in some way, so expectation values are always to be understood as averages over gauge configurations 
with some gauge-fixing term present. Which gauge is chosen makes no difference to any of the results in this paper, as all 
the improvement coefficients are gauge independent. 



(6) 




(7) 
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in terms of the fermion matrix M . Inserting the definition (||) into eq. we find 

S*{x,y) = [S{x,y) - ^A^(5(x - y)^ , (8) 

with improvement coefiicients 

b^^ = -\. (9) 
Av, = 1. (10) 

We use the lattice deha function 

5{x-y) = ^5^y, (11) 

where 5xy is the dimensionless Kronecker delta function. Although the propagators S and 5* in eq. (||) 
are of course gauge-dependent, the improvement coefficients 6^ and are gauge independent. 

Note that although we have used Kqw at intermediate stages in the discussion, our final result (||) 
only involves inverting the fermion matrix M . As noted before, Kqw is not well defined if Dq\y has 
eigenvalues exactly equal to 2/a, which will happen in topologically non-trivial configurations. On the 
other hand, the fermion matrix AI has no such problems. It is always invertible for mo > 0, so eq. (^ can 
be applied even in configurations with a non-trivial winding number. Our improved Green's functions 
will always be of such a form that only appears as a propagator. An identity equivalent to eq. (H) 

that will often prove useful is 

1 / o- ^ a 



Now that we have defined a propagator, we can look for an expression for the chiral condensate. The 
natural choice of order parameter for chiral symmetry is to take the trace of the improved propagator S^,, 

{ij{x)^(x))^ = TrS.{x,x)^^^Ty(^(l-^DGw)M-^) . (13) 

The final form of the improved chiral condensate is the same as that given in M . 



3 The fermion propagator in the free theory 

In this section we consider one explicit realisation of the Ginsparg- Wilson condition (P, namely Neu- 
berger's Dirac operator We show that the formulae we derived do indeed lead to results free of 0{a) 
effects. 

Starting from the massless Wilson fermion matrix Dy/i Neuberger introduces the matrix A, defined 

by 

A=l-aDw- (14) 
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It can then be shown that the operator 

DN=-(l-A/VAfA) (15) 



a 



satisfies the Ginsparg- Wilson condition (0) . 

In the free theory, the Wilson matrix is diagonal in momentum, and has the value 

Dwip) = i^+W, (16) 

where we use 

W = ^(1-cosap^), (17) 

^ = ^7^sinap^, (18) 

^sin^ap^. (19) 



s2 



Calculating Dn, the Ginsparg- Wilson matrix corresponding to Dw, we find Q 

A^A = (1 - + (20) 

and 

1 i^+(w-\ + V(i - + sA 

^^^P^ = ^ ,,,,, , -■ (21) 

If we expand Djsiip) for small p, we find 

I?w(p) =iy+^ap' + 0(aV), (22) 

so Dtv still has discretisation errors of 0{a). 

However, when we calculate Kn{p) according to the formula (^, we find 

Kn{p) = (23) 

al-W+ ^(l-I^)2+s2 

Finally in K^ip) we have an operator which anti-commutes with 75. When we expand it for small p, we 
find 

^ , . iE,7.(p.-ia^P^ + 0(aV)) 

^^(^^ = l-iaV + 0(aV) ^ ^ 
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Figure 2: The scalar part of the free propagator, (l/4a)Tr5(p), plotted against a^p^, for the various 
fermion matrices considered in this section. The Wilson propagator is the dotted line, the unimproved 
propagator from is the dashed line, the improved propagator Si, is the dot-dashed line, and the solid 
line is the continuum result too/(p^ + Wg). All curves are plotted for amo — 0.01. The momentum is 
taken in the direction (1,1,1,1). 

which has lattice errors of O(o^). 

However, does have problems of its own. It diverges when = and W > 1, which occurs at 
the 'doubler' momenta pa = (tt, 0, 0, 0), • • •. Because it has poles in momentum space, it would also be an 
extremely non-local matrix if Fourier transformed back into position space. This is, of course, inevitable: 
the Nielsen-Ninomiya theorem tells us that a fermion matrix that anti-commutes with 75 and has no 
doubling must be non-local. 

The improved propagator does not have these difficulties. It is 



S*{p, mo) 



^ + mo(l-W+ ^(1-1^)2 +52) /2 



2(w-l + ^y{l-W)^^ + 8"^^ /a2 +ml(l-W+ ^/{T^^WfT^'^ /2 
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= ^f±^ + 0(a^). (25) 

The denominator has no zeroes if mp > 0, and if mp = the only zero is at p = 0. The numerator 
vanishes at the 'doubler' momenta. If we expand this propagator as a power series in a, we see that there 
is no 0{a) term: the discretisation errors of 5* are of O(a^). 

In Fig. 1^ we compare the trace of the various lattice propagators with the correct continuum result. 
The Wilson propagator and the unimproved Ginsparg- Wilson propagator both deviate from the 

correct result when a?p^ « amo, but the improved propagator Si, is good up to a^p^ « 1. 



4 Improvement of flavour non-singlet fermion operators 

Next we want to improve the Green's function corresponding to a flavour non-singlet operator 

O = 4)Oil), (26) 

where O can include Dirac structure and covariant derivatives. Covariant derivatives can be represented 
by the usual formulae 

5^/ = ^^[U^{x)f{x + fi)-Ul{x-fi)!{x-fL)\, 

Id^ = ^[f{x + fi)Ul{x)-f{x-fL)U^{x~fi)], (27) 

or by any other expression with discretisation errors of 0{a^). We want our improved Green's function 
to be 

Gf^l/— O— ^ \. (28) 

a* \ Kgw + ma Kgw + mo / 

This will have discretisation errors of O(a^), because the propagators used have been improved, and the 
usual discretisation of O is valid up to terms of O(a^). 

We need a formula analogous to eq. (H), giving in terms of M~^. The simplest way to construct 
such an identity is by using eq. (^2|), which gives 



where 



O = {l + amQb^)(l~^DGw)o(l-^D. 



O + amob^O - ^(1 + amob^){DGwO + O Dgw) + -^{^ + amob^)DGwO Dgw- (30) 
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Equation (|29|) is written with the same wave function improvement factor as eq. because the quantity 
of physical interest is always the ratio G^/S"*, and we want the wave function factors to cancel. 



Equation ( p9|) is not the most general expression for . We can use the equation of motion 

1 



DgwM-^ = - ■ 



Too 



-M- 



1 — aTOo/2 ' 1 — aTOo/2 

to show that this is equivalent to the expression 

1 " 



Gy = 



1 + amob^ 



Go — 



where 



Go 
G° 



O + anioCoO — - 



O + amoCoO - -ci{DgwO + O Dgw) 



with 



Co 



Cl 



amoC2 



jamo 



1-ci 



4(1 , 
amoC2 



1 \7 ' 

-flTOo j 



1 — -^amo 2(1 — ^amoy 

1 — C2 — ^flTOo 

(1 - iaTOo)2 



(31) 



(32) 



-c2DgwO Dgw, 



(33) 
(34) 
(35) 



(36) 
(37) 
(38) 



There are two free parameters in this system of equations. The improvement coefficients ci and C2 can 
take any value, but once it is chosen, the values of the other improvement coefficients are fixed. This 
freedom comes from the equations of motion, which allow us to compensate for a change in one of the 
improvement coefficients by adjusting the other coefficients. For an example of this in the clover action 
see @. 

The terms proportional to G'^ and (O) can be interpreted as contact terms. The Green's function 
we are interested in has the form {^i'ipjOjkipk'4'i)- On the lattice we should expect to see a contact term 
of the form Sij{Ojki^k4'i) + {'4'i'4'jOjk)Ski, with a coefficient of 0(a), and a 'double contact term' of the 
form 5ij{0jk)Ski, with a coefficient of O(a^). 

If we are looking at the operator Green's function on-shell, these contact terms are irrelevant, because 
the fermion fields will all be well separated in position, and so the delta functions are all zero. In this 
case we can use any values of ci and C2, as long as we use the correct cq value. 

On the other hand, if we look at the off-shell Green's functions, we have to take the contact terms 
into account. One possibility would be that as well as computing Go one would also compute G^ and 
(O) and add them with the coefficients Xq and ijo- A more elegant procedure would be to use eq. (p9|). 



8 



i.e. to choose ci and C2 so that the contact terms are absent, and Go is improved both on- and off-sheU. 
The improvement coefficients for this special case are 



1 

ci = C2 = 1 - -amo, 

Co = -\- (39) 

Note that, as in the previous section, the matrix Kqw is only used heuristically. To find the improved 
Green's functions 5*^, and G^, the only matrix that we really have to invert is M . 

Finally, we mention another way to calculate the improved Green's functions S^, and of cqs. (^, 
(^). From eq. (^|) we can write the improved propagator as 

S. = ((i-^Z?gvv)m-i). (40) 

Since Dqw a-nd M commute, this can be rewritten in the more symmetric form 

= ( ( 1 - ^Daw) ' (l - ^Daw ) ■ (41) 

This formula gives exactly the same improved propagator as eq. (^, but it lends itself to a somewhat 
different interpretation. In eq. (H) we subtract an unwanted contact term present in the unimproved 
propagator, while eq. ( ^l|) can be interpreted as the propagator of a 'rotated' fermion field, which is 
similar to the picture of clover improvement presented in . 

Similar formulae can be written for the improved Green's functions, where instead of adding irrelevant 
higher dimension terms to O we perform a 'rotation' of the fermion fields. The resulting formula is 



(42) 



where 



O™* ^ (l - ^Daw) ' O (l - ^Daw) ' ■ (43) 

Using the identity (^2|) it is easy to see that eq. (|4^) is equivalent to eq. (p8|). The rotation approach 
requires us to take the square root of a matrix, which is not needed in the irrelevant operator approach, 
so it might be more costly to implement. 

As with the propagator, we can illustrate the effects of improving the Green's function by looking at a 
simple case in the free theory, namely the Green's function G^{p, nio) for the local vector current ■(/'7^tV'. 
Here t is a flavour matrix, Trr = 0, normalised so that Trr^ = 1. As in the case of the propagator, the 
most sensitive test is to look at the scalar part of the Green's function, which is shown in Fig. ||. We 
have chosen the quantity 



E 



i^^-i^Tr[G''(p,mo)r], (44) 
a 4o^ 



9 



10 



4^ Tr G'^T 



4a 



10 



10 



10 



10 



10 



: 1 1 

_ X 




: / 

y 

- 


X^"' unimproved GW 


- 


\ Wilson ■, \ ; 








\ ■ ^ - 
~\\ \ > : 




■\ ' - 
• \ 1 

\\ 1 
\ \ '.1 




\ \ '■. : 




improved GW ^, \ 




\ \ 
\ \ ■■. 




\ ; 




■^continuum- 




\ 



10" 



10 ^ 10 10 ^ 10 ^ 



10 



10^ lO'^ 



2 2 

a p 



Figure 3: The scalar part of the Green's function for the local vector current ip^/(4a^)Tr[G'^(p)T]. The 
Wilson fermion gives the dotted line, the result from Ginsparg- Wilson fermions with no improvement 
terms is the dashed line, the dot-dashed line is from the improved Green's function G^,, and the solid line 
is the continuum result 2tooP^/(p^ + tti^Y . As in Fig. || all curves are plotted for amo = 0.01, and the 
momentum is taken in the direction (1,1,1,1). 

which has the value 2tooP^/(p^ + Wg)^ in the continuum. Again, the Wilson action and the unimproved 
Ginsparg- Wilson Green's function deviate significantly from the desired continuum result at a^p^ « amo, 
while the improved Green's function only has errors of 0{a?) and remains reliable until a?p^ k, 1. 

As a final example in the free case, let us take an operator with a derivative, namely 



i - 1 
2 



(45) 



with fjL ^ V. This is an operator that can be used to compute the flavour non-singlet moment {x) of a 
hadronic structure function. In the continuum, the free Green's function for this operator has the value 



-iy -f mo 1 -ij 

7: (nt,Pu + iuPt,) T- 



mo 



(46) 
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Figure 4: The scalar part of the Green's function for the operator O^^ of eq. (|45|), i/(4a)Tr[G''''(p, mo)T]. 
The Wilson fermion gives the dotted line, the result from Ginsparg- Wilson fermions with no improvement 
terms is the dashed line, the dot-dashed line is from the improved Green's function G*, and the solid line 
is the continuum result 2mQp^p^ / {p^ + ™o)^- Fig.^ all curves are plotted for amo = 0.01 and the 

momentum is taken in the direction (1,1,1,1). 

The improved lattice Green's function found by applying the formulae in sect, ^is 



G'^''{p,mo) = S^{p,mo) — (7^sinap^ +7^sinap^) tS'*(p, mg). 



(47) 



where the explicit form of S'*(p, mo) is given in eq. (|25|). Since the discretisation errors of S** are 0{a^) 
it is clear that the difference between G'^^^^{p, mo) and Gt^ [p^ nio) is also of O(a^). This is illustrated in 
Fig. , where we plot the trace of G'"'. 
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5 Ward Identities 



In the continuum chiral symmetry relates the scalar part of the fermion propagator to the forward Green's 
function for the non-singlet pseudoscalar operator ip^^Tip: 

j5TS{p) + S{ph5T = 2mG^{p). (48) 

On the lattice we can check that the improved Green's fmictions of sect. || satisfy the same Ward identity: 

-f5TS,{p) + S^iph^T ^ 2moGUp) . (49) 

One might naively expect that this Ward identity would be violated at O(a^) by the discretisation errors, 
but in fact it holds exactly, even though both the left hand side and the right hand side have discretisation 
errors of O(a^). 

Note that all the Green's functions needed for Ward identities, like eq. (|49|), can be computed on 
the lattice. So we can hope to use them to determine the improvement coefficients in other improved 
theories, such as clover fermions, where we cannot calculate them analytically from first principles. 



6 Conclusions 

We see that using the Ginsparg- Wilson action there are many ways of constructing propagators and 
Green's functions which are correct to O(a^), both on-shell and off-shell. In this paper we have considered 
both, improvement by adding the Green's functions of irrelevant operators, and by rotating the fermion 
fields. 

The idea is very simple. In eqs. (|^) and (|2^) we have the Green's functions we wish to obtain, with 
the correct chiral symmetries, and therefore free from all 0(a) errors. However, there they are expressed 
in terms of the associated matrix Kq^^, defined in eq. (||). Fortunately, it is easy to find equivalent 
expressions in terms of the well-defined matrix M, so we never need to find or invert Kqw explicitly. 

The improvement coefficients are universal. It makes no difference which theory is being considered, 
Abelian or non-Abelian gauge theory, or a theory with no gauge symmetry, or what value of the coupling 
is being used. It does not even matter which operator is being considered. Operator improvement 
for other improved actions, such as clover fermions, will certainly not be so simple. The improvement 
coefficients in that case will certainly be functions of the coupling, and depend on the theory considered. 
Nevertheless, Ginsparg- Wilson fermions may give some hints as to what sort of improvement terms are 
needed, and what sort of contact terms are to be expected. 
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